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We present a numerical lattice quantum chromodynamics calculation of isovector form factors 
and the first few moments of the isovector structure functions of the nucleon. The calculation 
employs two degenerate dynamical flavors of domain- wall fermions, resulting in good control of chiral 
symmetry breaking. Non-perturbative renormalization of the relevant quark currents is performed 
where necessary. The DBW2 gauge action is used to further improve the chiral behavior while 
maintaining a reasonable physical lattice volume. The inverse lattice spacing, a -1 , is approximately 
1.7 GeV. Degenerate up and down dynamical quark masses of approximately 1, 3/4 and 1/2 times 
the strange quark mass are used. The physical volume of the lattice is about (1.9 fm) 3 . The ratio of 
the isovector vector to axial charges, gA/gv, tends to a somewhat lower value than the experimental 
value as the quark mass is reduced toward the physical point. Momentum-transfer dependences 
of the isovector vector, axial, induced tensor and induced pseudoscalar form factors are calculated. 
The Goldberger-Treiman relation holds at low momentum transfer and yields an estimation of the 
pion-nucleon coupling, g^NN = 15.5(1.4), where the quoted error is only statistical. We find that the 
flavor non-singlet quark momentum fraction (x) u -d and quark helicity fraction (x)Au-Ad overshoot 
their experimental values after linear chiral extrapolation. We discuss possible systematic errors 
for this discrepancy. An estimate for transversity, {l)$ n -5d = 0.93(6) in MS at 2 GeV is obtained 
and a twist-3 polarized moment, di, appears small, suggesting that the Wandzura-Wilczek relation 
holds approximately. We discuss in detail the systematic errors in the calculation, with particular 
attention paid to finite volume, excited state contamination, and chiral extrapolations. 

PACS numbers: 11.15.Ha, 11.30.Rd, 12.38.Aw, 12.38.-t 12.38.Gc 



I. INTRODUCTION 

In this paper we report numerical lattice quantum 
chromodynamics (QCD) calculations of the isovector 
form factors and moments of the isovector structure 
functions of the nucleon on gauge field configurations 
with two degenerate flavors of dynamical domain-wall 
fermions (DWF) 0,0, Q generated by the RIKEN-BNL- 
Columbia (RBC) Collaboration. 

Four isovector form factors parameterize neutron (3 de- 
cay: the vector and induced tensor form factors from the 
vector current, 

(P\^(x)\n) = u p [i^gviq 2 ) - q\<y\^gT{q 2 )} u n e iq - x , 

and the axial and induced pseudoscalar form factors from 
the axial current, 

(p\ A t( x )\ n ) = Up bv75ffA(<7 2 ) + iq^j59p{q 2 )] u n e lq ' x , 

(2) 
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which are given here in the Euclidean metric conven- 
tion. Thus, q 2 as used in this paper stands for Euclidean 
four-momentum squared, and corresponds to the time- 
like momentum squared as q 2 M — — q 2 < in Minkowski 
space. We also note a slight change in the sign con- 
vention of the pseudoscalar form factor from our ear- 
lier publications, such as Ref. Q. Here q — p n — p p 
is the momentum transfer between the proton (p) and 
neutron (n). In the limit \q\ — > 0, the momentum trans- 
fer should be small because the mass difference of the 
neutron and proton is only about 1.3 MeV. This makes 
the limit q 2 — ► 0, where the vector and axial form fac- 
tors dominate, a good approximation. Their values in 
this limit are called the vector and axial charges of the 
nucleon: gv = gv(q 2 = 0) and gA = gA{q 2 = 0). Exper- 
imentally, gv = cos 9c (with the Cabibbo mixing angle 
6 C ), and g A = 1.2695(29) x g v @. 

These form factors are calculable on the lattice [H, 0, 
i, ©, IHHHI but quite often are prohibitively 
complicated if one uses conventional staggered or Wilson 
fermions. The staggered fermions, with their complicated 
flavor/taste structure, make even definition of baryon op- 
erators difficult. The Wilson fermions make the necessary 
current renormalization complicated due to large explicit 
violation of chiral symmetry. The domain-wall fermions 
(DWF) [T3 |. with their exponential suppression 

of chiral symmetry breaking, make such renormalizations 
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easy. In particular, currents that are connected by chi- 
ral transformation such as vector and axial ones should 
share the same renormalization. Thus, the ratio Qa/ Qy 
is naturally renormalized in DWF lattice calculations H . 

Earlier (mostly quenched) lattice calculations [HI, [l6|, 
\v?\ of this ratio gave significant underestimations by up 
to 20%. In a quenched calculation with DWF we found 
much of this deficit comes from the small lattice volumes 
employed in those earlier studies [1]. With the suffi- 
ciently large lattice volume of about (2.4 fm) 3 used in 
our quenched calculation with DWF and the improved 
DBW2 gauge action [ID, [l9[ , we observed that the ratio 
does not depend strongly on the quark mass and ob- 
tained a value of 1.212 ± 0.027(stat) ±0.024(norm) in the 
chiral limit. The value is almost consistent with experi- 
ment, and combined with the very weak dependence on 
the quark mass, does not require any chiral extrapolation 
but a linear one in quark mass. Thus an obvious question 
arises: do these observations hold when the quarks are 
treated dynamically? 

There has been recent interest in the pseudoscalar form 
factor, gp(q 2 ), due to the MuCap Collaboration's new 
high-precision experiment studying ordinary muon cap- 
ture (OMC) by protons, /jTp -> v^n ;2fj]. The OMC ex- 
periments determine the induced pseudoscalar coupling 
gp = m fl gp(q 2 ) at q 2 = 0.88m 2 , where m M is the muon 
mass. The new experiment yields gp = 7.3±1.1, which is 
consistent with the value predicted by chiral perturbation 
theory, g£ hPT = 8.26 ± 0.16 [HI], but much smaller than 
the last world average for OMC, (s£ MC )avc = 10.5 ± 1.8 
given in Ref. [12] , and the value from a single experiment 
of radiative muon capture (RMC), fi~p — > v^ny, g P MC = 
12.4±1.0 [13]. Recently, a quenched DWF calculation @ 
reported a result, gp — 8.15 ± 0.54(stat) ± 0.16(norm), 
consistent with the MuCap experiment. It is interest- 
ing to see what value is obtained in a rif — 2 dynamical 
DWF calculation. 



In this paper we also report on this form factor, gp(q 2 ). 
It is the induced part of Eq. |J2J|, directly related to 
the pion-nucleon coupling, g^NN, which should satisfy 
niNgA = Fn9nNN) a t finite momentum transfer, 



2m N g A {q 2 ) - q 2 gp{q 2 ) 



q 2 + m 2 



(3) 



with a more traditional convention of = ~ 
92 MeV. These should hold up to residual dependence on 
the momentum transfer q 2 ; the pion-nucleon coupling at 
high momentum transfer may significantly differ from its 
value near zero momentum transfer. On the lattice there 
may also arise some corrections from the finite lattice 
volume. 

Since we study the dependence of these form factors at 
relatively low momentum transfer, we can try to extract 
corresponding mean square radii, (r )v,t,a, defined by 

gv,T,A{q 2 )/gv,T,A{Q) ~ 1 - \q 2 {r 2 )v,T,A H • 
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The structure functions are measured in deep inelastic 
scattering of leptons from nucleons 0, [H, [27], [H, 
H M M, HI IM H HI, the cross section of which is 
factorized in terms of leptonic and hadronic tensors: 



da 



2M 



a 



"°d 3 fc' s-M 2 (q 2 ) 2 
The leptonic tensor is known to be 



lfj,u{k,k ) — 2 ^k^k u -\- h^h v g^ 



(4) 



(5) 



Hence, the cross section provides us with structure infor- 
mation about the target nucleon through the hadronic 
tensor, W^ Vl which is decomposed into symmetric unpo- 
larized and antisymmetric polarized parts: 



J 



W^(x,Q 2 ) = ( - 



q^q h 



Fl {x,Q 2 )+(p»-^ 



q- 



F 2 (x,Q 2 ) 



(6) 



W^(x,Q 2 ) = ie^q p (^( gi (x,Q 2 ) + g 2 (x,Q 2 )) 



q ■ SP a 2 

2 — 92{X, Q ) 



V = q- P, S 2 = -Af 2 , and x = Q 2 /2v, and Q 2 = \q 2 \. 
The unpolarized structure functions are F\(x, Q 2 ) and 
F 2 (x,Q 2 ), and the polarized, gi(x,Q 2 ) and g2(x,Q 2 ). 



Their moments are described in terms of Wilson's op- 
erator product expansion: 
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2 f dxx^F^x^Q 2 ) 
Jo 

[ dxx n ~ 2 F 2 (x,Q 2 ) 
Jo 



2 / dxx n g x {x,Q 2 ) 
Jo 



2 / dxx n g 2 {x 1 Q 2 ) 
Jo 



E cg( M 2 /Q 2 )5 ( M )) (x n )M+0(l/Q 2 ), 

q—u^d 

E 41(M 2 /Q 2 ,.9(M))<x") g (^) + 0(l/Q 2 ), 

f=u,d 

E e &(/* 2 /Q 2 ,<7M) 0Oa 9 M + o(i/g 2 ), 



1 n 



2n 



j E [4«(^ 2 /Q 2 ^(^))^(M)-2<„(M 2 /Q 2 ,.9(M))(^>A 9 (/i)] +o(i/g 2 ), (8) 



q—u,d 



r 



where ci, c 2 , &x, and e2 are perturbatively known Wilson 
coefficients and (x n ) q (fi), (x n )^ q (fj,) and d n ([i) are calcu- 
lable on the lattice as forward nucleon matrix elements 
of certain local operators. 

Again, the conventional staggered or Wilson fermions 
would complicate such lattice calculations for the same 
reasons as discussed for the form factors. The DWF cal- 
culations are simpler because of easier renormalizations 
due to good chiral symmetry. In particular the first mo- 
ments (x) u —d (quark momentum fraction) and (x) Au-Ad 
(quark helicity fraction) share a common renormaliza- 
tion, and so their ratio is naturally renormalizcd in DWF 
calculations [3^ |. 

During the last few years lattice computations have 
provided many interestin g r esults for these structure 
function moments [4j, M, l a , 1 1 71 . l36l . 37 1 38l 39 

0, H EE EE S3, fiafia M M m M m 




both quenched and full QCD. These calculations provide 
first-principles values for the moments of structure func- 
tions at leading twist. One of the major unresolved issues 
in these previous calculations is significant overestima- 
tion of the moments {x) u — d and (x) Au—Ad compared with 
results from fits to the experimental data [171, [36j, |38|, |48| . 
In our quenched calculation [36| of these quantities we 
encountered similar overestimations. Interestingly, how- 
ever, the ratio of these quantities, (x) u —d/ {%) Au—Ad 
showed very weak dependence on the quark mass and 
agreed well with the experimental ratio in the chiral limit. 
Again an obvious question is whether this behavior sur- 
vives when the quarks are treated dynamically. 

We address these questions with a lattice QCD calcu- 
lation using two degenerate dynamical flavors of domain- 
wall fermions and DBW2 rectangle-improved gauge ac- 
tion [57[. Three bare sea quark masses of 0.02, 0.03 and 
0.04 in lattice units (corresponding to pion masses of 
about 0.5, 0.6 and 0.7 GeV) are used with about 200 
gauge configurations each. The lattice cutoff is a -1 ~ 
1.7 GeV, and the spatial volume is about (1.9 fm) 3 . The 
lattice cutoff is sufficiently high to allow us to take advan- 
tage of various benefits of DWF such as good chiral and 
flavor symmetries in performing fully-nonperturbative 
renormalization for our nucleon observables. Unlike our 



past quenched calculations 4, 36] which did not show sig- 
nificant dependence on quark mass in most of the observ- 
ables, we sometimes observe significant deviations at the 
lightest quark mass from the heavier quark mass results. 
The pion mass times lattice extent, m^L, is slightly less 
than 5 at the lightest quark mass, which may be problem- 
atic for calculations involving a large hadron such as the 
nucleon. We must therefore perform quark-mass extrap- 
olation on a case-by-case basis and with great caution. 

The rest of the paper is organized as follows. In Sec- 
tion [Til we nrs t briefly summarize the numerical ensem- 
bles used in the calculation. Then we discuss the choice of 
nucleon source and sink and the operators used to calcu- 
late the form factors and moments of the structure func- 
tions. Finally we briefly summarize the numerical non- 
perturbative renormalization methods. In Section [Ml we 
begin with discussion of the form factors at zero momen- 
tum transfer. We show the conservation of the isovec- 
tor vector current is under good control and then give 
an evaluation of the nucleon axial charge. Then we dis- 
cuss momentum dependence of various form factors. In 
particular the Dirac mean-squared charge radius is ex- 
tracted. The Goldberger-Treiman relation is shown to 
hold and provides an estimate of the pion-nucleon cou- 
pling. Finally, we describe our calculations of some low 
moments of the structure functions. We summarize our 
conclusions and future plans in Section [TVl 



II. NUMERICAL METHOD 

We use DWF lattice ensembles generated by the RBC 
Collaboration These were generated with two de- 

generate flavors of dynamical quarks described by the 
domain-wall fermion [l], 0, [H, 03 action and gluons de- 
scribed b y th e doubly-blocked Wilson (DBW2) gauge ac- 
tion [3[igl. There 

are three such ensembles with sea 
quark masses of m soa = 0.02, 0.03, and 0.04 in lattice 
units, which respectively correspond to about 1/2, 3/4 
and 1 times the strange quark mass, or degenerate pseu- 
doscalar meson masses of about 500, 600 and 700 MeV. 
The inverse lattice spacing is about 1.7 GeV, set from the 
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p-meson mass, yielding a physical volume of (1.9 fin) 3 . 

We study nucleon matrix elements using quark prop- 
agators obtained with Gaussian gauge-invariant sources 
with radius (r gauss ) of 4.35 for all three quark masses and 
additionally with Gaussian smearing radius of 8 for the 
lightest quark mass. The former parameter was previ- 
ously used for the neutron dipole moments calculation 
on the m soa = 0.03 and 0.04 ensembles [58]. The lat- 
ter is a better size for overlap with the ground state 
of the nucleoli. For details of our source definition see 
Ref. [HI]. We use only unitary valence quark mass values, 
"Valence = ^sea = wi/ j when computing quark propaga- 
tors. We use two different sequential sources generated 
with source-sink locations, [£ S rc , ^sink] = [0,10] and [15, 
25], where the source-sink separation i sep is fixed as 10 
lattice units, on a given gauge configuration for all three 



quark masses with the smaller 



As for the larger 



Rgauss i we use a different choice of source-sink separa- 
tion, t scp = 12. The longer separation causes the larger 
statistical noise. Therefore, we use three sources placed 
at [tsrc^sink] = [0, 12], [9, 21] and [18, 30] to increase the 
statistics. Later, we will discuss the possibility of excited- 
state contaminations in our calculations by comparing re- 
sults obtained from two sets of parameters, {t scp , r gauss }, 
on the lightest quark ensemble. The parameter values 
and the number of configurations used from each ensem- 
ble are summarized in Table |TJ We also compile some 
basic physics results from Ref. 57| in the same table. 

We define our nucleon two-point functions with a nu- 
cleon interpolation field x an d smearing parameters A 
and B as follows: 



(XA(t,m B (0, -P)) = YPoufo s)<P, s\x B \0)e- Eip)t +■ = mY *l' f+mN VzA(p)z B ( P ) e"^' + • • ■ (9) 



r 



with normalized states defined as where E(p) = yj m? N + p 2 . The ellipsis in Eq. ((9J denotes 

excited-state contributions, which can be ignored in the 



(0\xa\p,s) = y/z A (p)u s (p), (10) case oft >1. 

and spinors satisfying 1 + 7t 

We often use the projector V + = — - — - to the 

u s (p)u s (p) = E{p)^f t — 17 • p + mjv, (11) positive-parity states, 



r^WsrcSp) = \^(V + (xA(t s ^,P)XB(t SI c^p))) = m ^ E mN VZA(P)Z B (P) e - g ^(«.c fc -W) + . . .. (12) 



These two-point functions provide appropriate normal- the nucleon mass, 
ization factors when we extract matrix elements from 

three-point correlation functions, as well as estimates for Now we define the three-point functions: 



G M ,AB(tsrc,i,tsnk) = (XA (*snk, Psnk) J ji (t , q) Xb ( f src , Psrc)) 



+ irf • f + m N _ E>f-i1J +mN e -*(w-*) e -*C«») + . . . ( 13 ) 



2E' M 2E 



where q = f - p, p = p SIC , f? = p snk , E = E(p STC ) and 
E' = -E(psnk)- The operator is appropriately selected 
for each observable of interest discussed in the following 
subsections. 



We calculate the nucleon isovector form factors and 
some low moments of structure functions, namely the 
isovector vector (gv), induced tensor (gr), isovector ax- 
ial (<7a); and induced pseudoscalar (gp), form factors at 
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TABLE I: Summary of simulation parameters and the numbers of configurations. The pion mass, the nucleon mass and the 
pion decay constant (F n ) are also tabulated. 



m f 


0.02 


0.03 


0.04 


^snk ^src 


12 10 


10 


10 


^src 


0,9,18 0,15 


0,15 


0,15 


^"gauss 


8 4.35 


4.35 


4.35 


# of conf. 


185 each 220 each 


220 each 


220 each 


(GeV)[57j 


0.493(5) 


0.607(4) 


0.695(4) 


m N (GeV)[57] 


1.28(2) 


1.43(2) 


1.55(2) 


Fjv (GcV)f57] 


0.1141(7) 


0.1232(7) 


0.1329(7) 



both zero and finite momentum transfer, and structure 
function moments corresponding to the momentum frac- 
tion (x) u -d, helicity fraction (x)au-a<i, a twist-3 moment 
di, and transversity (l)su-Sd- The respective choice of 
the lattice operators for these observables are the same 
as in earlier RBC reports 0, H, [H, [59[ on quenched cal- 
culations, and are briefly summarized in what follows. 

Before we move on to discuss the operator choices 
for our observables of interest, we would like to briefly 
explain our different choice of source-sink separation, 
i scp , and Gaussian smearing radius, r ga u SS , for the light- 
est quark mass ensemble (see [H, |4£| [fiO] for the ear- 
lier calculations). We started our calculation by setting 
^scp = 10 time units and the Gaussian source smear- 
ing radius to 4.35, parameters previously used for the 
neutron dipole moments calculation on the m/ = 0.03 
and 0.04 ensembles 58]. Although the new results of the 
Gaussian source in the heaver quark-mass region agree 
with the heavier quark-mass points of the box source, for 
the same calculation with m/ = 0.02, we find a signifi- 
cant decrease in the axial charge compared to our previ- 
ous box-source calculation, which used a separation of 12 
time units as shown in Fig.[TJ This discrepancy prompted 
us to examine whether there is a systematic error due to 
excited-state contamination in the new results. 

The effective masses plotted in Fig. [5] show that the 
Gaussian source with smearing radius 4.35 has a ground 
state plateau that begins at t — t SIC = 7 or 8 (top figure) , 
while the box source appears to plateau sooner (bottom 
figure). Thus one may suspect that excited-state con- 
tamination accounts for the effect observed in Fig. [T] 
Since the major difference between the two calculations 
is the source type and separation, we extend the source- 
sink separation to t sop — 12 for the Gaussian source (the 
same as in the previous box-source calculation) and also 
increase its radius to 8 in order to reduce excited-state 
contamination. Indeed, the Gaussian source with radius 
8 provides similar quality of plateau in the effective mass 
plot, as shown in the middle panel of Fig. [2j However, 
as shown in Fig. [TJ the resulting g^l Qv with the second 
parameter set {i S cp, ?~gauss} = {12,8} is consistent with 
that of the first one {10,4.35} rather than the previous 



result of the box source. So, the discrepancy between the 
results is probably not caused by excited-state contami- 
nation. 

As we mentioned before, to reduce statistical fluctu- 
ations, we use multiple sources in the present studies. 
This is another difference from the previous calculation, 
where only a single source was utilized. To make this 
point clear, we show the dependence of bare gv and gA 
on the location of the current insertion t for each choice 
of the source location in Fig. [3J In the top panels, all 
calculations of the bare gv show clear plateaus between 
source and sink locations. Although the slight depen- 
dence on the source location may be observed among the 
three sources in the {t scp , r gauss } = {12,8} case, all the 
values of these plateaus agree with one another within 
the statistical errors. We see no evidence of excited-state 
contamination in the bare gv- 

On the other hand, for the bare gA (bottom panels), 
both the results obtained from box and Gaussian sources 
with t scp — 12 exhibit larger fluctuations and less clear 
plateaus, while two consistent plateaus clearly appear in 
the cases of t scp = 10. This indicates that the larger 
separation causes the larger fluctuations, since the abso- 
lute values of three-point functions are exponentially sup- 
pressed as a function of t sep . Indeed, the t scp = 12 Gaus- 
sian source results obtained from a single source-sink lo- 
cation [isrc^sink] = [9,21] agree with those of the box 
source due to the large statistical errors, while the aver- 
age from all three sources for {^sep 1 Rgauss 

} = {12, 8} pro- 
vides good agreement with the results of {t S cp 7 Rgauss} = 
{10,4.35}. 

In this context, the lower value of the bare gA ob- 
tained from the present studies is more statistically sig- 
nificant, since the box source results have been obtained 
only for a single source. Although a systematic er- 
ror from excited-state contaminations in the result of 
{^scp, Rgauss} = {10,4.35} might be hidden, owing to the 
large fluctuations in that of {t scp , r gauss } = {12,8}, we 
may choose to live with the former for the final result. 
This is mainly because the finite momentum calculation 
suffers much from such large fluctuations. We obtain re- 
sults for the form factors only from the case of t sep = 10. 
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in terms of the dimensionless vector Gv and induced 



FIG. 
with 



m; [GeV ] 

1: Renormalized isovector axial charge g_\/ . 
box [431 and Gaussian smeared sources 



jv calculated 
with various 



source-sink separations (t scp ). The value from the Gaussian 
source with t scp = 10 and smearing radius 4.35 at the light- 
est quark mass of m/ = 0.02 deviates significantly from the 
previous box-source calculation. In that case the box size is 
chosen for better overlap with the ground state of the nucleon 
and the longer source-sink separation (i sop = 12) is adopted 
to avoid possible excited-state contaminations. This discrep- 
ancy is not resolved by using a Gaussian source of an radius 
8 and the same source-sink separation. As discussed in the 
text, this is likely caused by the larger statistical fluctuation 
in the longer source-sink separation. 



tensor Gt form factors. Here the E M in Eq. (fl3|) is 
l»G v tf)-q^^. 

Likewise the isovector axial current in Eq. (QJ is rewrit- 
ten as 



(p\A+(0)\n) = u p 



G P (q 2 ) 



with the dimensionless axial Ga and induced- 
pseudoscalar Gp form factors. Note that the latter 
is normalized with twice the nucleon mass 2m^ unlike 
in some muon-capture literature where it is normalized 
with the muon mass m M . This is for numerical conve- 
nience in this paper, where the nucleon mass calculated 
on the lattice is heavier than its physical value. Here 
the E„ in Eq. (SJ is W G A (q 2 ) + iq^ 5 ^l. (In the 
above, gv,A(q 2 ) = G v ,A{q 2 )-) 

The right-hand sides of the above two equations have 
the most general form consistent with Lorentz covariance. 
The momentum transfer q — p n — p p becomes very small 
in the forward limit, because of the small mass difference 
between the neutron and proton. 

We use two projection operators to help us extract the 
momentum dependence of the form factor: 



A. Form factors 

For numerical convenience, we redefine the form factors 
as dimensionless quantities. Thus, the isovector vector 
current in Eq. ((!]) is rewritten as 



<p|V+(0)|n) = % 



l^G v (q 2 ) - <?aCTa m 



GAq 2 ] 
2mjy 



(14) 



-Tt{VG^AB (*snk, t] p src , Psnk)), 



^|t,As(^ snk ' ^PsrcPsnk) — T^*"' v '/' 

(16) 

where we choose V = Va, = ~P+j5jz for both axial and 
vector currents and V = V+ for vector current. We com- 
bine three- and two-point functions to remove time de- 
pendence and redundant source or sink normalization (z) 
factors 14 ill: 



p __ ,GG(*sro t, f 8 nk)Psrc)Psnk) / T^q (t, t sn k, Psrc ) Tgg (i S rc , t, Psnk)^LG^ STC > ^ snk ' Psnk) \ ,.. _> 

j " ~ r (2) i, I z s x 1 r (2) 7 , z x r (2) , z ^ F (2) j ' z r I ( - i '> 



where L denotes a local (point) source or sink and G 
stands for a Gaussian-smeared one. In this work, we 
fix the sink momentum to zero. Therefore, in the axial- 
current case, Eq. (fT7|) gives 



n Al = 



i 



^2E(p)(E{p) + m N ) 
S iz (E(p) + m N )G A ■ 



PzPi 

2mjy 



Gi 



(18) 



with i £ {x, y, z} for corresponding insertion of operator 



Ai. In the vector current case, we use the projection 
operator P + on V4 and Pa z on V\ and V% to extract Gy 
and induced Gt] the necessary equations are respectively 



' E(p) + m N ( q _ E(p) - m N ^ 



2E(p) 



-i Py 



^/2E(p)(E(p) + m N ) 

+iPx 

y/2E(p)(E(p) + m N ) 



2th 
(Gy + Gt), 

(G v + Gt). 



(19) 
(20) 
(21) 
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FIG. 2: Nucleon effective mass plots from two-point functions. Gaussian source with smearing radius, r gauss = 4.35 (left), 
(center), and box source (right) from the calculation in Ref. [43| . 
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FIG. 3: Comparison of the bare isovector vector charge gry (top) and axial charge gA (bottom) from source-sink separation 
i SC p = 10 Gaussian source (left panels), t scp — 12 Gaussian source (middle panels) and i scp = 12 box source (right panels). 



One can in principle solve for Ga{q 2 ), Gp(q 2 ), Gv(q 2 ), 
and Gt(<? 2 ) with overconstrained data for each q 2 . 

B. Moments of structure functions 

Of the moments of the structure functions summarized 
in Eq. (j5J), we calculate those which do not require finite 
momentum transfer: the quark momentum fraction (x) q , 
helicity fraction {x)& q , transversity (l)gq and the twist-3 
moment d% . Of these we calculate the isovector contribu- 
tion for the former three: q = u — d. This simplifies their 
renormalization as the flavor-singlet contribution drops 
out. Thus we fully nonperturbatively renormalize these 
three so the results are compared with the experiments 
under an assumption of good isospin symmetry. On the 
other hand we do not renormalize the twist-3 moment: 
here our interest is to look at its individual up and down 
components to see if the Wandzura-Wilczek relation (6lT | 
holds. Our choice of operators (E M in Eq. (fiU)) ) for the 
three-point functions follow an earlier RBC paper (3o| 



and are summarized in Table HT1 



C. Nonperturbative renormalization 

In order to compare our calculation with the experi- 
mental values, we need to establish the proper connec- 
tion to the continuum through renormalization. For- 
tunately, the well-preserved chiral and flavor symme- 
tries of the domain-wall fermions (DWF) make this task 
much easier than in the cases of more conventional 
fermions [H, [H, [H . 

For the form factors, the chiral symmetry of DWF as- 
sures the isovector vector and axial local currents which 
are used in the present calculation share a common renor- 
malization: Zy = Za, up to higher order discretization 
errors, 0{a 2 ). Further, since the vector current is con- 
served, the vector renormalization is calculated from the 
vector charge, gv, as Zy = 1/gv- These will be demon- 
strated in the numerical results section. Note the induced 
tensor and pseudoscalar form factors share this common 
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TABLE II: Operators used in our structure-function mo- 
ments calculations, including the notation for the operator, 
the explicit operator form, the hypercubic group representa- 
tion, the correlator ratios and the projection operators used 
in the numerical nonperturbative renormalization in Eq. (|26[) . 



quark momentum fraction {a;)^ 



O q i4 = q 
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hypercubic group rep. 
correlator ratio 



a 
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Pt 



= m N {x) q 



NPR projection V\ A 1 = j 4 p 4 - | J2i=i,3 

quark helicity fraction (x)& q 



hypercubic group rep. 
correlator ratio R (x)i 



7» D 4 +74 D-i 
67 



r o q 

°3pt 



NPR projection 



= m N {x)A q 

; «75 (73P4 + 74P3) 



transversity {l)s q 



°ll = 9750-34(? 



hypercubic group rep. Q{ 

] 
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R Wi q 



correlator ratio 
NPR projection 



r o aq 
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°3pt 



c> 



pt 
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NPR projection 



73 Da -74 D-i 



Rd 



r 0' jq 

n '"[34] 
°3pt 

C2pt 



^[34] = ^75 (73P4 ~ 74P3) 



renormalization as well and do not require any additional 
calculation for renormalization. 

For structure function moments, we follow Ref. [66| 
and implement the regularization-independent 
momentum-subtraction (RI/MOM) scheme to non- 
perturbatively renormalize the operators. This was first 
applied to structure functions by Gimenez et al. [67l ]. 
The procedure may be described as follows. First, 
the Fourier transform of the Green function with our 
operator of interest is constructed: 



G 0r (p; a) = ^2 ' 



-ip-(x-y) 



<V(*)Or(<Wl/)>, (22) 



x.y 



where Or is one of 02 4 , OH, or °[34] for the first- 

moment operator renormalization. In the lattice calcu- 
lation, we need to produce Fourier-transformed point- 
source S (pa; 0) and point-split-source D fl S(pa; 0) propa- 



gators, 



S(pa;0) = J2e- ipx S(x;0) (23) 

X 

D^S(pa;0) = £ le*** [S(x; -ftU^-fi) 

X 

- S(x;fi)Ul(0)], (24) 



to construct the Green function in Eq. (|22|) . In this work, 
we use momenta (p x ,p y ,p z ,p t ) ranging from (0, 0, 0, 0) to 
(2,2,3,5) in units of 2-k/L, 

The next step is to truncate the external legs from the 
Green function: 

A 0r (p;a) = S(p;a)- x Go r (p;p';a)S(p';a)-\ (25) 

where Z RI is obtained after the projection 

Z 0r (n;a,y 1 Z q (p;a) = ^Tr (A 0r (jj; o)flr)| p 2 =/1 i. (26) 

The projection operators for various structure function 
operators are listed in Table HU Note that p must fall 
inside the renormalization window Aqcd <C p -C l/a. 

The following steps allow us to convert the renor- 
malization constants into the continuum MS scheme 
at 2 GeV. See Appendix Q] for specific details on the 
strong-coupling scheme-matching coefficients, RGI run- 
ning equation and anomalous dimensions for various op- 
erators. 

1. Obtain Z Kl (p): 

The ratio of Zo r (p;a)/Z q (/j,;a) to ZA/Z q (p,;a) is 
computed and yields Zo r (p; cl)/Za- Each of the 
factors in the ratio is first exptrapolated to the chi- 
ral limit, nif — — m rcs , at fixed momentum. Us- 
ing Zy = l/.9v, we can determine Zo T (p;a) in 
Eq. (f2"6"]l . the renormalization constant in the RI 
scheme, which we denote as Z Kl . (See the lightly- 
filled circles in Fig.|4j) 

2. Convert to MS scheme: 

We are interested in continuum quantities, mostly 
calculated in the MS scheme. The conversion 
factors between RI and MS schemes for the op- 
erators discussed here have been calculated in 
Refs. [H, [H. To get Z MS (p), we use the three- 
loop a s (p) running coupling constants defined in 
Ref. [6| with A(n/ = 2) = 300 MeV0. 

3. Remove (ap) 2 lattice artifacts: 

To remove 0((ap) 2 ) errors that might mimic con- 
tinuum scale dependence, we first divide out the 
continuum running factor. The resulting renormal- 
ization factors should be scale- independent, so any 
remaining (ap) 2 dependence must be a lattice arti- 
fact. Hence, Z MS ' KGl can be obtained from a fit to 
the form / = A(pa) 2 + B. (See Fig. gj) 
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FIG. 4: Renormalization constants in the chiral limit. The 
lightly-filled circles are the renormalization constants in RI- 
MOM scheme, and the diamonds are MS-scheme at 2 GeV. 
The fits (solid lines) are used to remove residual (ap) 2 arti- 
facts. 



4. Running to 2 GeV: 

Finally, we use the RGI formula to obtain 
Z m (2 GeV). 

Table Mil summarizes the renormalization factors in MS 
scheme at 2 GeV for each operator. These renormaliza- 



TABLE III: Summary of renormalization factors in MS 
scheme at 2 GeV in the chiral limit. The number of con- 
figurations is about 100. 



— TOrcs 1.28(7) 1.23(7) 0.776(14) 



tion constants will be applied in Section lill CI 



III. NUMERICAL RESULTS 
A. Vector and axial charge 

The isovector vector and axial charges, gv and gA, are 
defined as the zero-momentum-transfer limits of the cor- 
responding form factors, Gy(q 2 ) and Ga(<? 2 )- Because of 
chiral symmetry, a Takahashi-Ward identity ensures that 
the two currents, which are related by chiral transforma- 
tion, share a common renormalization: Za = Zy up to 
a lattice discretization error of 0(a 2 ,rnj a 2 ). Since the 
vector current is conserved, its renormalization is eas- 
ily obtained as the inverse of the vector charge gv = 
Gv(,Q 2 — 0). Thus, by calculating the ratio of the three- 
point functions for gA/gvi we get the renormalized axial 
charge, (g A y en \M- 

Let us first discuss the vector charge. We note, how- 
ever, that since the vector charge is conserved, it can not 
provide a transition matrix element from one state to an- 
other, in particular from any excited state to the ground 
state. Thus the precocious plateau does not necessarily 
mean that excited state contamination is absent in the 
signal after one time unit. As mentioned earlier, to avoid 
excited-state contributions at the lightest quark mass, 
we have lengthened the source-sink separation to 12 time 
units and used a larger smearing radius for the Gaussian- 
smeared source and sink. From the fitting ranges of 3-9 
for t sep = 12 and 3-7 for i scp = 10 (see also Fig. [5]), we ob- 
tain the plot in Fig. [51 This figure shows that larger fluc- 
tuations appear in the case of the longer time separation. 
By linearly extrapolating three points of t sep = 10 to 
zero quark mass, we obtain an estimate of gv — 1.308(9) 
and Zy = 0.764(5). The resulting value agrees reason- 
ably well with the value of Za = 0.75734(55) found in 
Ref. [57| . which was determined from the ratio of two- 
point meson correlation functions. Thus we conclude 
that excited state contamination in our value of gv is 
negligible. 

In Fig. [7] we present the raw results for the bare ax- 
ial charge. After taking an average of all sources, the 
a plateau emerges, even in the case of t scp — 12 (top 
figure), albeit with larger statistical errors. As can be 
seen from this figure, the plateaus settle in after about 
three time slices for both t sep cases at the lightest quark 
mass. This indicates that excited-state contamination is 
negligible compared to the statistical error. Correspond- 
ing renormalized values are obtained from the ratio of 
correlation functions, and are presented in Fig. [8] As ex- 
pected, we also find that plateaus settle in after about 
three time slices for all cases. By using the fitting ranges 
shown in this figure, we obtain the values in Table IIVI 
and Fig. H 

We observe that our axial charge value at the lightest 
quark mass deviates significantly from the heavier quark 
mass points. This deviation causes the linear extrapo- 
lation to the physical point to deviate from experiment 
by more than five standard deviations. (The obtained 
value is given as gA/gv — 0.89(6) with a poor value of 
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FIG. 5: Bare vector charge as a function of the current in- 
sertion time. Results are averaged over all sources. The 
top figure depicts t scp — 12 calculation at the lightest quark 
mass, rrif = 0.02. The other three figures are obtained from 
iscp = 10 calculations at all three quark masses (in order of in- 
creasing mass from top to bottom). The fit ranges are shown 
by horizontal lines. 
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FIG. 6: The isovector vector-charge (<7v) bar ° as a function of 
the pion mass squared. 



There is a caveat that we have to arbitrarily fix at least 
one of the four parameters in SSE. However, if one con- 
siders only the finite-volume correction SgA (which is 
gA(L) — 9a{L = oo)) within this model, SgA depends 
only on two phenomenological parameters, the AT- A cou- 
pling (ca) and A-A coupling (gi) , which can be barely 
fixed. Furthermore, if the iV-A coupling is set to zero, 
the Beane- Savage formula for the finite- volume correc- 
tion to gA reduces to that of HBChPT at leading order. 
In the original paper [7l|, SgA is predicted to be positive. 
This remains true in the case of ca = 0, where SSE re- 
duces to LO HBChPT. In this sense, this formula fails to 
account for the negative SgA observed in the quenched 
DWF calculation. 

In Ref. [zl] , a finite- volume study of gA has been done 
with rij = 2 dynamical clover simulations. They also ob- 
serve negative SgA in both quenched and dynamical sim- 
ulations. Although their lattice data points are outside 
of the range of applicability of SSE (or LO HBChPT), 
they can fit their data using the SSE formula for gA with 



X 2 /dof ~ 11.3) Omitting the lightest mass point, a con- 
strained linear fit gives gA/gv — 1.23(12) at the physical 
point, in good agreement with experiment. As mentioned 
111 Section HH our physical volume is about (1.9 fin.) . In 
light of our previous quenched calculation [4[ where we 
saw that a small volume resulted in a small axial charge, 
one may worry that a similar finite-volume effect occurs 
here. 

Since our quenched calculation [4| was published, a way 
to interpret the finite-volume effect has been proposed 
by Beane and Sava ge J7ll| within the small scale expan- 
sion (SSE) scheme [72|], which is one possible extension 
of HBChPT with explicit A degrees of freedom. The 
SSE scheme can provide milder quark-mass dependence 
for the axial charge, which seems to be consistent with 
lattice results, while either leading order (LO) or next- 
to-leading order (NLO) HBChPT shows strong quark- 
mass dependence in the vicinity of the chiral limit [73| . 



TABLE IV: Bare isovector vector and axial charges, g\? Te 
and an d their ratio (su/gvO barc which is by-definition 

renormalized, extracted from data summarized in Figs. [S] [7] 
and [8] 



m f 


"{^sep j ^"gauss } 


bare 

9v 


_barc /_ / \rcn 

9 a (gA/gv) 


0.02 


{12, 8} 


1.28(3) 


1.28(12) 1.01(10) 


0.02 


{10, 4.35} 


1.299(6) 


1.27(4) 0.97(4) 


0.03 


{10, 4.35} 


1.291(3) 


1.53(3) 1.19(3) 


0.04 


{10, 4.35} 


1.2865(13) 1.51(2) 1.171(16) 



1 The correspondence between the SSE notation and the param- 
eters used in Ref. |7l| can be found in Ref. [74|. Here, we prefer 
to use the original SSE notation. 
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FIG. 7: The bare axial charge as in Fig. [5] 



FIG. 8: The axial-to- vector charge ratio as in Fig. [5] 



the finite-volume corrections, which should be negative. 
They insisted that the parameters adopted in the Beane- 
Savage paper do not fully reflect decoupling constraints, 
which are guaranteed to reduce SSE to LO HBChPT 
when Cj\ — 0. Indeed, their adopted parameter set yields 
5g A < 0. 

Apart from the concept of SSE, this suggests that 
the sign of SgA predicted from SSE is very sensitive 
to adopted parameters. We find that the choice of c A 
strongly affects the sign of SgA, as shown in Fig. [Til] We 
fix the chiral limit value of the axial charge gA, the pion 
decay constant F„ and the N-A mass splitting An to 
their physical values and set g\ = 3.0 as in Ref [71j |. 
We also set L = 1.9 fm. We vary the remaining pa- 
rameter of the iV-A coupling ca from to 1.5. Here, 
we recall that ca = 1.5 is adopted in Ref. [74|, while 
ca = 1 in Ref. [7l| . The predicted correction at our sim- 
ulation point, less than 5%, is negligible. Furthermore, 
the sign of SgA seems to be positive, against our expec- 
tation. Therefore, the finite-volume correction from SSE 
can not account for our data either qualitatively or quan- 
titatively. Again, it is likely that the heavy quark masses 
in our simulation do not allow use of such formulae. We 



will come back to this finite-volume question when we 
discuss the momentum dependence of the form factors 
and quark momentum and helicity fractions. 



B. Momentum dependence of form factors 



As described in the previous subsection, data sets 



with {t sep , 



3 } = {10,4.35} do not seem to suffer 



from excited-state contamination. Therefore, we focus 
on those data sets for analysis of the form factors at fi- 
nite momentum transfer in this subsection. 



1. Vector current 

Under the present assumption of exact isospin symme- 
try that ignores strange quark contribution, a straight- 
forward exercise in SU(2) Lie algebra leads to the follow- 
ing well known relations between the weak vector-current 
matrix elements in neutron beta decay and the isovector 
part of the nucleon electromagnetic current matrix ele- 
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where N represents p (proton) or n (neutron). The 
isovector ones G" E M are defined by G V E M = G P E M — 

r~in 

As such we can compare our lattice results computed at 
finite momentum transfer with the experiments. The eas- 
iest quantity for such a comparison is the charge radius, 
which is conventionally defined as (r 2 ) = —QdFi/dq 2 . 
The proton mean-squared charge radius is known ex- 
perimentally Q as 0.7656(119) fm 2 and the neutron, - 
0.1161(22), resulting in an estimate of 0.636(12) fm 2 for 
the mean-squared "Dirac" radius of the nucleon. (For 
details, see Appendix A of Ref @). 

The raw results from the lattice are presented in 
Fig. [TTJ From these we extract the values of the Dirac 



FIG. 9: The axial-to- vector charge ratio as a function of the 
pion mass squared. The lightest mass points with either t scp = 
10 or t scp = 12 deviate significantly from the heavier mass 
points. Omitting the lightest mass point, a constrained linear 
fit is performed. 




FIG. 10: Finite volume correction of the axial charge <?a is 
predicted within 0(e 3 ) SSE. The sign of SgA is very sensitive 
to the choice of ca- 



ments: 

{p\u la d\n) = (p\u-f a u - d la d\ P ) = (j}\f a m \p) ~ (H-O) 

(27) 

where j® m = ^u-f a u — ^djad. This relates the weak 
vector and induced tensor form factors with the isovector 
part of the nucleon electromagnetic form factors: 



F? (q 2 ) - F?tf) = G v (q 2 ) 
F%{q 2 )-FZ{q 2 ) = G T (q 2 ). 



(28) 
(29) 



These are respectively are called the Dirac and Pauli form 
factors. They are related to the more conventional Sachs 
electric Ge and magnetic Gm form factors: 



N, 



GZtf) - W) + W) 



(30) 
(31) 
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FIG. 11: Bare Dirac isovector form factor as a function of the 
current insertion time at the heaviest pion mass (m/ = 0.04). 
The horizontal lines represent the average value (solid lines) 
and their one standard deviations (dashed lines) in the fitted 
range. 

isovector form factor, normalized by the value at q 2 = 0, 
as listed in Table |V] In Fig. Q21 we plot all three 
raf results for Gv{q 2 ) as a function of Euclidean four- 
momentum squared together with the dipole form of 
(1 + q 2 /AI 2 )~ 2 with the empirical value of the isovector 
Dirac dipole mass M — 0.857(8) GeV, which is evalu- 
ated from the electric charge and magnetization radii of 
the proton and neutron. Different symbols represent the 
values obtained from different quark masses. As can be 
seen, there is no large quark-mass dependence. Our cal- 
culation points are located far from the empirical curve. 

The slope of the form factor at q 2 — determines the 
mean-squared radius, which can be related to the corre- 
sponding dipole mass M as (r 2 ) = 12/M 2 . To extract 
the mean-squared radius from our data, we simply adopt 
the dipole form for fitting three lower q 2 points including 
the q 2 = value. The fitted values of the Dirac dipole 
mass My and corresponding mean-squared Dirac radii 
{r v ) are listed in Table ED 

We show the quark-mass dependence of measured 
mean square Dirac radius in Fig. 1131 Although it is dif- 
ficult to extrapolate the lattice estimates at the three 



TABLE V: Momentum transfer dependence of the isovector 
form factors, n 2 refers to the momenta carried by the vertex 
in units of (2ty/L) 2 . 



m f n 2 = 


n 2 = l 


n 2 = 2 


n 2 = 3 


n 2 = 4 


0.02 1.299(6) 0.86(3) 0.68(3) 0.63(7) 
G v 0.03 1.291(3) 0.882(18) 0.73(3) 0.62(4) 
0.04 1.2865(13) 0.930(12) 0.725(17) 0.61(3) 


0.51(17) 

0.48(5) 

0.58(5) 


0.02 N/A 
G T 0.03 N/A 
0.04 N/A 


2.13(15) 
2.67(13) 
2.81(9) 


1.53(11) 
2.05(10) 
2.04(7) 


1.28(18) 1.1(4) 
1.57(12) 1.11(13) 
1.69(8) 1.18(13) 


0.02 1.27(4) 
G A 0.03 1.53(3) 
0.04 1.51(2) 


1.09(4) 
1.16(2) 
1.21(2) 


0.92(5) 
0.97(4) 
0.99(3) 


0.87(9) 
0.87(5) 
0.88(4) 


0.58(19) 

0.69(8) 

0.82(8) 


0.02 N/A 
G P 0.03 N/A 
0.04 N/A 


9.3(8) 
9.5(6) 
9.7(5) 


4.9(5) 
6.1(4) 
6.1(4) 


3.3(5) 
4.2(4) 
4.4(3) 


1.8(7) 
2.6(4) 
3.1(4) 
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FIG. 12: Momentum-transfer dependence of vector form fac- 
tors. A dashed curve corresponds to the dipole form with the 
empirical value of the dipole mass My = 0.857(8) GeV, which 
is evaluated from the electric charge and magnetization radii 
of the proton and neutron. 



heavy pion masses to the physical pion mass, the mild 
mf dependence in our observed range of m 2 allows us 
to perform a simple linear extrapolation. Then we ob- 
tain the mean-squared radius of the vector form factor at 
the physical point: (ry) = 0.281(32) fm 2 , which is much 
smaller than the experimental value of 0.636(12) fm 2 . 

However, it is well known that this particular quan- 
tity, as well as the pion charge radius, has a logarith- 
mic divergence in heavy baryon chiral perturbation the- 
ory (HBChPT) as one approaches the chiral limit (75|. 
The dotted curve plotted in Fig. Q21 is the expected 
quark-mass dependence of (r v ) in the vicinity of the 
physical point within leading order (LO) HBChPT. (See 
Appendix 1 2 bl for more details.) As can be seen, the 
very steep m 2 dependence, which is associated with the 
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TABLE VI: Dirac mean-squared charge radius obtained from 
fitting to the conventional dipole form. 



m f 


M v (GeV) (H>)(fm 2 ) 


0.02 


1.40(5) 


0.239(19) 


0.03 


1.49(4) 


0.209(12) 


0.04 


1.57(3) 


0.190(8) 


Expt. 


0.857(8) 


0.636(12) 
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FIG. 13: The mean-squared Dirac radii (r v ) from Gv(q 2 )- 
A simple linear extrapolation with respect to the pion mass 
squared underestimates the experimental point (asterisk). A 
dotted curve is the predicted quark-mass dependence of mean- 
squared Dirac radius from LO HBChPT in the vicinity of the 
physical point. 



logarithmic divergence, near the physical point is pre- 
dicted by HBChPT at leading order. This may account 
for the smaller values of our measured mean-squared ra- 
dius, which are calculated in the heavy quark mass region 
(m, > 0.49 GeV). 

Thus, we may expect that there is non- linear behav- 
ior in terms of m 2 in the vicinity of the physical point. 
This expectation may look like somewhat contradictory 
to what we argue for gA with HBChPT and SSE. How- 
ever, we recall that the logarithmic divergence can not 
be cured by higher-order corrections, whereas the strong 
quark-mass dependence of gA, which simultaneously im- 
plies slower convergence of the chiral expansion, should 
be modified by higher-order corrections. Indeed, both 
the one-loop effective field theory approach with explicit 
A degrees of freedom (SSE) [z3 and two loop HBChPT 
calculation [76j present flat quark-mass dependence at 
least down to ~ 0.3 GeV. To observe this, we have to 
at least push the pion mass down to 0.3 GeV, which is ap- 
proximately the location of the intersection between the 
LO HBChPT curve and the extrapolated line obtained 
from our three data points in Fig. Q21 We leave this issue 
to future calculations. 

Next, we extract the induced-tensor (Pauli) form fac- 
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tor, Gt, from the vector current; the raw data is listed 
in Table fVl and an example of nif — 0.04 at various 
transfer momenta is shown in Fig. [14] This corresponds 




FIG. 14: The bare Pauli isovector form factor as in Fig. [TT] 



to the combination Ff — F% of the nucleon electromag- 
netic form factors. The renormalized Pauli form factor 
can be obtained from G^ n {q 2 ) = G T (q 2 )/G v {0) with 
vector current renormalization Zy = 1/ Gv (0) . Fig. 1151 
shows a comparison of our data and the dipole form 
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FIG. 15: Momentum-transfer dependence of the induced- 
tensor form factors. The dashed curve corresponds to the 
dipole form with the experimental value of the dipole mass 
Mt = 0.778(23) GeV and nucleon magnetic moment F 2 (0) = 
3.70589. 

G5? n (0)/(1 + q 2 /M T ) 2 with the experimental values of 
G5? n (0) = F 2 (0) = 3.70589 and M T = 0.778(23) GeV, 
as described in Appendix A of Ref @ . The heaviest two 
pion mass points are almost degenerate, as was the case 
with the Dirac form factor. However, the data from light- 
est pion mass here is closer to the experimental values. 
Again, this could either be interpreted as a trend to- 
ward the experimental values with decreasing pion mass 
or merely finite- volume effects. 

Furthermore, we fit our G^ n (q 2 ) to the dipole form to 



extract the Pauli mean squared radius, which is related to 
the corresponding dipole mass Mt by (r T ) = 12/M 2 ,. In 
contrast to the dipole fit on Gv(q 2 ), it is a two-parameter 
fit since we do not have data on the value of G!j? n (q 2 ) at 
q 2 = without the q 2 extrapolation. Here, the value 
of the G!j? n (0) is associated with the difference between 
the proton and neutron magnetic moments y, v — [i n = 
1+F2(0). The fitted results are summarized in Table IVHl 
Here we perform the dipole fit with the two lower q 2 
points. 



TABLE VII: Pauli mean-squared charge radius obtained from 
fitting to the conventional dipole form and the extrapolated 



value of the ratio G v M (q 2 ) /G v E {q 2 ) 
q 2 fitting form. 



at q 



with the linear 





F 2 (0) Mr (GeV) (r^fm 2 ) 


GV(0)/G|(0) 


0.02 
0.03 
0.04 


2.57(37) 1.29(16) 0.28(7) 
2.89(22) 1.53(12) 0.20(3) 
3.32(17) 1.36(6) 0.25(2) 


3.69(41) 
4.23(26) 
4.24(17) 


Expt. 


3.70589 0.778(23) 0.773(32) 


4.70589 



In Fig. \W\ we plot values of (r 2 -) and F^O) as a func- 
tions of pion mass squared along with a naive linear 
extrapolation to the physical pion mass. In the upper 
panel of Fig. |T6l we linearly extrapolate through all three 
points, finding (r 2 -) = 0.17(11) fm 2 . This is much smaller 
than the experimental measurement, 0.773(32) fm 2 , and 
similar to what we observed for the Dirac mean-squared 
radius. Lighter pions in future measurements would be 
desirable to see whether there is an increase in radius as 
we approach the physical pion mass. 

In the lower panel of Fig. [Till although the value mea- 
sured at the heaviest point is close to the experimental 
value, the extrapolated value of F^O) at the physical 
point also tends to somewhat underestimate the experi- 
mental value. However, we recall that the value of F2(0) 
is highly dependent on our adopted fitting form. For 
example, a monopole fit yields a larger F2(0) which lies 
closer to the experimental value. Therefore, our estima- 
tion of F 2 (0) should carry a large systematic uncertainty 
due to the q 2 extrapolation. 

We also have an alternative way to evaluate F 2 (0) = 
fi p — ji n — 1. The ratio of the isovector magnetic form fac- 
tor G v M {q 2 ) and the isovector electric form factor G E (q 2 ) 
provides the difference between the proton and neutron 
magnetic moments in the forward limit, fj, p — fi n = 
G v M (0)/G v E (0), which is related to 1 + F 2 (0). Experi- 
mentally, it is known that this ratio shows no q 2 depen- 
dence at low q 2 since both form factors are well fitted 
by the dipole form with the comparable dipole masses. 
Therefore, this ratio may have milder q 2 dependence than 
G T T n (q 2 ) @]. Some RBC results on related ratios for the 
heavier quark masses were reported in Ref. [58| . 

In Fig. [T71 we show the q 2 dependence of the ratio 
G\, I (q 2 )/G E (q 2 ), which clearly exhibits mild dependence. 
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FIG. 16: The mean-squared Pauli radius (r^) (top) and the 
value of G^ n (0) = ^2(0), which is associated with the dif- 
ference between the proton and neutron magnetic moments, 
— fJ-n = 1 + ^2(0), (bottom) and as a function of pion 
mass squared. A simple linear extrapolation with respect to 
the pion mass squared underestimates the experimental point 
(asterisk). 



We therefore use a simple linear fitting form with respect 
to q 2 for an alternative evaluation of the value /i p — fx n . 
In Fig. [TS1 we plot values of [i v — /i„, from two deter- 
minations, 1 + F 2 (0) and Gm(0)/Gs(0), as functions of 
pion mass squared. Both determinations are roughly 
consistent with each other within the statistical error. 
Weak m 2 dependence is observed; the heavier two points 
in the case of Gm(0)/Ge(0) are the same within error, 
while slight downward dependence appears in the re- 
sults of 1 + i*2(0). The simple linear extrapolation of 
Gm(0)/Ge(0) yields \i v — fi n = 3.4(7), which gives the 
better agreement with the experimental value. We may 
quote this value for our final value of \i p — since the 
latter approach appears to have smaller systematic error 
in the ^-extrapolation. 



2. Axial-vector current 

In Fig. [BJ] we present the raw data for the axial form 
factor at finite momentum transfer. These allow us to ex- 




qlGeV 



FIG. 17: The q 2 dependence of the ratio G v M (q 2 )/G v E (q 2 ). A 
simple linear extrapolation with respect to q 2 is utilized for 
an alternative evaluation of fip — fi n thanks to its mild q 
dependence. 
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FIG. 18: Comparison of two determinations of fi p — fi n . 
There is no appreciable m^. dependence for the results from 
G V M (0) / G V E (0) . A simple linear extrapolation is applied to 
them. The resulting extrapolated value at the physical point 
shows the better agreement with the experimental value (as- 
terisk) . 



tract the momentum-transfer dependence of the form fac- 
tor as listed in Table W\ and plotted in the top of Fig. [501 
The axial form factor is phenomenologically fitted with 
the dipole form, at least at low q 2 ; so are the Dirac and 
Pauli form factors 2l| . 



The dashed curve in the top 
of Fig. [201 shows the dipole form with an experimental 
value of the axial dipole mass M A = 1.026(21) GeV 
Our lattice data lie above the experimental curve. Notice 
that the likely finite volume effect observed earlier for the 
axial charge leads to the non-monotonic behavior of the 
form factor with quark mass, and leads to the downward 
curvature of the mean-square axial radius. 

To extract the mean-squared radius of the axial vector 
form factor, we perform the dipole fit to the three lowest 
q 2 points and convert the axial dipole mass Ma into the 
mean-squared axial radius (r\) 



The obtained values are 
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FIG. 19: The bare axial form factor as in Fig. Qj] 



TABLE VIII: Dipole fits for momentum-transfer dependence 
of the axial form factor. 



m f Ma (GeV) (r\) (frn 2 ) 

0.02 2.26(25) 0.091(20) 

0.03 1.72(7) 0.158(11) 

0.04 1.90(7) 0.129(9) 



Expt. 1.026(21) 0.444(19) 



listed in Table ECCD and the bottom of Fig. M We also 
plot the mean-squared axial radius as a function of the 
pion mass squared. As can be seen, the value of (r\) 
at the lightest pion mass shows a large reduction beyond 
statistical fluctuations. This large reduction in the axial 
radius suggests finite-volume effects could be significant 
at the lightest quark mass. 

There is renewed interest in the pseudoscalar form fac- 
tor, Gp(q 2 ), due to the recent MuCap Collaboration [2(3] 
high-precision experiment studying ordinary muon cap- 
ture by protons, fi~p — > u^n, and because of improved 
electroweak radiative correction calculations [77( that al- 
low precise extraction of the form factor from these ex- 
periments. Using the new MuCap results, a value of 
g P = 2^G P (0.88m2) = 7 . 3±L1 is found |o|,[27l which, 
using p(5ac and chiral perturbation theory, is now in 
good agreement with QCD, g P = 8.26 ±0.12 [Hj]. Ear- 
lier, a TRIUMF group obtained a value of 12.4± 1.0 
Clark et al. [78[ got gp — 10.6 ± 1.1 after reanalyzing the 
TRIUMF data. Including the new MuCap result, the 
"world average" is 8.7 ± 1.0 [z3]. 

The induced pseudoscalar form factor is obtained 
from Eq. (fig]) , and is shown in Fig. [JU This allows us 
to extract its momentum dependence which is listed in 
Table [V] and plotted in Fig. [22] As is noted in the pre- 
vious section, we use a dimensionless definition for this, 
as well as the other three form factors, and use twice 
the nucleon mass 2m n estimated for each ensemble of 
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FIG. 20: (Top) Momentum-transfer dependence of vector 
form factors. The dashed curve corresponds to the dipole 
form with the experimental value of the dipole mass Ma = 
1.026(21) GeV, evaluated from pion electroproduction exper- 
iments. 

(Bottom) Mean-square axial radii (r A ) from Ga(q 2 )- The 
shaded regions are naive linear extrapolations (with errors) 
of our data points to the physical pion mass. An unexpected 
reduction of the measured {r\) at the lightest quark mass is 
observed. 



configurations. This is different from some experimental 
analyses where the muon mass is used for normalization. 
Momentum-transfer dependence of the pseudoscalar form 
factor Gp(q 2 ) has so far only been studied by one pion 
electroproduction experiment at Saclay in 1993 [79| with 
range below 0.2 GeV 2 . 

Fig. [JU shows the q 2 dependence of the ratio 
of the induced pseudoscalar to axial form factors, 
Gp(q 2 )/GA(q 2 )- Though our lattice momenta transfer 
are much higher than experimentally explored values, 
they trend upward toward the experiments. 

Another topic regarding the induced pseudoscalar form 
factor is its relation with the axial form factor through 
the PCAC relation and pion-pole dominance (PPD) 
model Hi: 



G P P V) 



Am 2 N G A (q 2 



(32) 



q* + 

To see how the pion-pole behavior is preserved in mea- 
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FIG. 21: The bare pseudo-scalar form factor as in Fig. 1111 
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FIG. 22: Momentum-transfer dependence of the induced 
pseudoscalar form factor compared with experiments. 



m=0.02 

■ m^O.04 

+ m=0.06 

3|£ experiment 



_l I I I I I I I I I I I I I I I I I L 



q [GeVl 

FIG. 23: The ratio of G P {q 2 ) and G P PD (q 2 ) as a function 
of four-momentum squared q 2 . Deviation of the ratio from 
unity indicated by dashed line is the deviation from the PPD 
model. However, the deviation is rather modest. At least, the 
q 2 dependence of the PPD model approximately accounts for 
that of measured Gp(q 2 ). 
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sured Gp(q 2 ), we consider the following ratio: 

G P (g 2 ) 
Gp PD (g 2 )' 

which is inspired by the above PCAC prediction. If the 
measured Gp(q 2 ) has exactly the same form described in 
Eq. (j32|) , this ratio yields unity in the entire q 2 region. 

In Fig. [531 we plot the above ratio as a function of four- 
momentum squared q 2 . There is no appreciable q 2 depen- 
dence. Thus, the q 2 dependence of the PPD model ap- 
proximately accounts for that of measured Gp{q 2 ). Four 
different q 2 points of the ratio cv PPD reveal a q 2 inde- 
pendent plateau within the statistical errors. We simply 
take the weighted average of cv PPD within four measured 
q 2 points, then plot them against the pion mass squared. 
As shown in Fig. [Ml the average values of a PPD gradu- 
ally approach unity as the pion mass decreases. A simple 
linear extrapolation of the quenching factor a PPD to the 
chiral limit yields 0.85(8), which is only about 2 stan- 
dard deviations away from the theoretical expectation. 



FIG. 24: The quenching factor a PPD is plotted as a function 
of the pion mass squared. 



One can obtain gp with input from a PPD : 

_ 9A 2mjyTO M 

9P a ^(l + q 2 /M 2 A ) 2 ( m l + q 2 y [M > 

where q 2 = 0.88m 2 with m M = 0.10568 GeV and — 
0.13957 GeV. M A and g A are 1.38(16) GeV and 1.23(12) 
respectively from naive linear fit using two heavier pion 
mass points. We found 

g P = 7.68 ± 1.03 (35) 

with statistical error evaluated by jackknife analysis. 
This is consistent with the "world average" value of 
8.7 ±1.0 [z3. 

The pseudoscalar form factor Gp is related to the 
pion-nucleon coupling, g^NN, through the Goldberger- 



18 



Treiman relation [8111: 



2m N G A (q 2 ) - q' 



G P (q 2 ) _ 2g vNN F 7r ml 
2m n q 2 + ml 



(36) 



At zero momentum transfer, the relation reduces to the 
following form: 



m N9A = F-xg-nNN- 



(37) 



Testing the relation experimentally in this zero- 
momentum form is a little tricky because the value of 
the pion-nucleon coupling, g^NN, is not known at zero 
momentum transfer, but at the pion pole. Worse, the 
pion-pole value of the coupling varies from analysis to 
analysis. A conventional partial- wave analysis of the 
pion-nucleon elastic scattering data to 2 GeV provides 
13. 75 ±0.1 5 [12] while a more recent reanalysis in Ref. [83j 
gives g-nNN — 13.169 ±0.057. Substituting the latter, we 
find there is a slight discrepancy from the Goldberger- 
Treiman relation at zero momentum transfer, Aqt, [84j | 
defined as 



A G t = 1 - rn N g A /F^g v 



NN 



(38) 



which is (2.259 ± 0.591)% if we use the mean nucleon 
mass m N = 938.9 MeV. 

Testing the relation at finite momenta transfer is more 
difficult, as the pion-nucleon coupling is even more poorly 
known from experiment. Thus it is useful to extract the 
pion-nucleon coupling, g^NN, from our lattice data: 



9nNN,liLt 



^AT,lat5A,lat 



lat 



(39) 



See Fig.[25]for g^NNM* results. A linear fit is performed 
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FIG. 25: Lattice estimate for the pion-nucleon coupling, 
SvrJVJVi from the Goldberger- Treiman relation as defined in 
Eq. (|39[) . Omitting the lightest mass point, a constrained 
linear fit is performed. 

without the lightest mass point to give an estimate of 
9ttNN — 15.5 ± 1.4 at the physical point. This should be 
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FIG. 26: Momentum-transfer dependence of the Goldberger- 
Treiman discrepancy as defined in Eq. (I40|) . 



compared with the experimental estimates at the pion 
pole such as 13.75 ± 0.15 or 13.169 ± 0.057 in the above. 

Now let us examine the Goldberger- Treiman relation 
by looking at the discrepancy Agt(<? 2 ) defined by the 
lattice quantities as follows: 



1 - 



2m,N 



lat 



257rATAf,lat^7r,lat'7l^ lat 



(40) 



see Fig. for the value of A GT (q 2 ). Note again that we 
are using the pion-nucleon coupling value determined at 
zero-momentum transfer, and cannot account for its vari- 
ation at finite momentum transfer. From our data for the 
two heavier quark mass values, we observe the discrep- 
ancy is strongly dependent on the momentum transfer 
beyond about q 2 k, 0.5 GeV 2 . This suggests the pion- 
nucleon coupling at such high momentum transfer values 
is very much different from the low momentum-transfer 
region. Note our pion mass squared is roughly m 2 , = 
0.25, 0.36 and 0.49 GeV 2 , or below the lowest q 2 in the 
plot, and the physical pion mass is even lower. Thus our 
data also suggest the Goldberger- Treiman relation holds 
at low momentum transfer, q 2 < 0.5 GeV 2 . In that re- 
gion the lightest quark mass result is in broad agreement 
with the heavier masses, albeit with large statistical er- 
ror. It is desirable to investigate this further at lower 
momentum transfer. 



C. Momentum and helicity fractions 

Let us now turn our attention from the elastic form 
factors to the deep inelastic structure functions. As is 
well known and has been summarized in the introduction, 
what are calculable on the lattice in regard to the struc- 
ture functions are their low-order moments. In this work 
we limit ourselves to the lowest order non-trivial mo- 
ments that are calculable with zero momentum transfer. 
As in previous RBC reports, the operator renormaliza- 
tions at the chiral limit are obtained nonperturbatively 
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using the RI/MOM scheme as described in Sec. Ill Cl 

We first discuss the quark momentum fraction, (x) u -d, 
or the first moment of the unpolarized structure func- 
tion. The bare lattice three-point correlators are shown 
in Fig. [27| and listed in Table [fx] The results for t scp = 12 
and t scp = 10 at the lightest pion mass differ by about 
two standard deviations. This may suggest excited state 
contamination in this quantity, even though its plateau 
looks very nice and flat over the range 3 < t < 7 for 
^scp = 10. From the top panel in Fig. [23 the plateau 
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FIG. 27: Bare momentum fraction as a function of the op- 
erator insertion time. Results are averaged over all sources. 
The top figure depicts the t scp — 12 calculation at the lightest 
quark mass, m/ = 0.02. The other three figures are obtained 
from t sep = 10 calculations at all three quark masses (in order 
of increasing mass from top to bottom). The fit ranges are 
shown by horizontal lines. 
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FIG. 28: Renormalized first moment of the unpolarized struc- 
ture functions or the quark momentum fraction, (x) u -d- Note 
that the circles indicate the data at each pion mass with two 
errorbars: the inner are statistical errorbars and the outer are 
errorbars propagated with renormalization factors. The solid 
square is the extrapolated point; the star is the experimen- 
tal value; and the band shows the chiral extrapolation with 
jackknife-calculated uncertainty. 



for t sop — 12 is not so clear or flat, owing to the much 
larger statistical fluctuations in the data for this larger 
source-sink separation. However, since it is low on the 
left and high on the right, the average value is insensitive 
to the choice of (a symmetric) fit range. In the earlier cal- 
culation, using a box source on a single time slice and a 
source-sink separation of 12 [43|,[6(|, the momentum frac- 
tion was consistent, within large statistical errors, with 
the t sep = 10 result. In this study, we can not conclude 
whether this is a true systematic effect or just statistics, 
especially given the quality of the i scp = 12 plateau. 

We show the results for all the measurements but only 
use the result from t sep = 10 at the lightest pion mass 
for extrapolation. Similar results would be obtained us- 
ing the i scp = 12 result at the lightest mass since the 
larger error on this point does not tightly constrain the 
fit. Using the fit ranges indicated in the figures and the 
renormalizations presented in subsection III C[ we arrive 
at the values of (x) u -d shown in Fig.[28j Note that there 
is only small pion-mass dependence; therefore, we use a 
simple linear fit. The quality of this fit in terms of x 2 
per degree of freedom is reasonable at 0.28. The chiral 
extrapolated value is (x) u -d = 0.282(19), overshooting 
the experimental value, 0.154(3) [13, [H, [3f|, by five to 
six standard deviations. Here, the trend of the i sep = 12 
data point is suggestive, but again, not conclusive be- 
cause of the relatively large statistical error. 

The quark helicity fraction, (%)Au-Adi is the- first mo- 
ment of the polarized structure function. The bare three- 
point correlator values are shown in Fig.[2H]and the values 
are listed in Table IIXI Similar remarks and conclusions 
about the plateaus for < scp = 10 and 12 as for the mo- 
mentum fraction hold here as well. The renormalization 
factors at the chiral limit for the helicity fraction can be 
found in Sec. Ill Cl Fig. I3H1 shows the renormalized quan- 
tities at each pion mass point. A linear fit to these values 
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FIG. 29: The bare helicity fraction; symbols as in Fig. 1271 



versus gives a chiral extrapolation of 0.286(25), over- 
shooting the experimental value of 0.196(4) l2<| HI] 
by almost four standard deviations. The quality of the 
fit is poor, as x 2 per degree of freedom is 6.78. In con- 
trast to the momentum fraction, the result at the light 
mass deviates from the heavier pion mass points, even at 
t sep = 10, a situation similar to the deviation observed in 
the axial charge. This trend is welcome since the experi- 
mental value apparently must be approached from above 
and suggests trying fits different than linear, to which we 
turn below, after discussing the ratio of the momentum 
fraction to the helicity fraction. 

In an earlier quenched calculation [36} we reported 
that the ratio of the momentum and helicity fractions, 
(x)u—d/{%)Au—Ad, agrees well with the experimental 
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FIG. 30: Renormalized first moment of the polarized struc- 
ture functions or quark helicity fraction, (iJau-aj, as a func- 
tion of m\. (Similar graphics description as in Fig. 1281 ) 
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FIG. 31: Renormalized ratio {x)Au-Ad/ {x) u -d of the quark 
momentum and helicity fraction. A linear fit is shown with an 
error estimate at the physical pion mass. Experiment gives 
0.78(2). 



value, 0.78(2). The quark mass dependence of this ratio 
is presented in Fig. [3U Unlike the quenched case where 
the ratio is only weakly dependent on the quark mass and 
is in rough agreement with experiment, there is a notice- 
able deviation for the lightest pion mass. A naive linear 
extrapolation gives very poor fit quality, % 2 /dof = 15.7, 
and a value at the physical pion mass of 0.93(7), which 
is about two standard deviations above the experimental 
one. This may be due to systematic effects like finite vol- 
ume, or excited state contamination; the t sep = 12 point 
is consistent with experiment and shows the weak mass 
dependence observed in the quenched case. 

Since our pion masses are rather heavy (700, 600 and 
500 MeV), it is not surprising that the chiral extrapo- 
lation to the physical pion mass (140 MeV) sometimes 
misses the experimental values. Also our heavy pion 
masses probably invalidate the application of a leading- 
order chiral form. The chiral extrapolation should ulti- 
mately be studied with more realistic lattice QCD ensem- 
bles such as the 2+1-flavor ones being generated jointly 
by RBC and UKQCD Collaborations!! H|, a work 
which is in progress. Nevertheless, it is a worthwhile 
exercise to try to extend the extrapolations beyond lin- 
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TABLE IX: Bare quark momentum and helicity fractions and 
their naturally renormalized ratio at each sea quark mass. 



m f 


tscp 




{x) u - d /{x) A u-Ad 


0.02 


12 


0.195(17) 0.22(3) 


0.99(14) 


0.02 


10 


0.236(9) 0.255(11) 


0.93(4) 


0.03 


10 


0.242(6) 0.306(6) 


0.790(17) 


0.04 


10 


0.253(4) 0.306(6) 


0.827(12) 
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FIG. 32: Renormalized quark momentum fraction, (x) u -d, 
and helicity fraction, (x)Au-Ad- The phenomenological fit 
(brown band) is according to Refs. [87, S&, \S§^ . The dashed 
lines are the linear extrapolation for comparison. (Graphics 
conventions as in Fig. 1281 ) 



ear forms. We now discuss how some such attempts fair 
with our data. 

We begin with a parametrization by Chen, et al. [13, 
HH Hil (the relevant equations are summarized in ap- 
pendix r2~cl ) Their prescription is to fix three parameters, 
the axial charge ffA,oxp, pion decay constant -F-^exp and 
the scale [i (set to the pion mass), to their experimental 
values and fit the coefficients C and e(/i 2 ). The results 
are summarized in Fig. 1321 The fit yields momentum 
fraction 0.147(9) after extrapolation, which is consistent 
with the experimental value. However, the quality of 
the fit is poor, % 2 /dof = 2.19; this is caused by a strong 
downward curvature in the fit form due to log terms. The 
corresponding helicity fraction extrapolation results in a 
value of 0.170(14), again, consistent with the experimen- 
tal value. The fit for the helicity fraction, while still poor, 
is improved over the linear one, with a \ 2 P er degree of 
freedom of 3.57. 



Another alternative is to consider the dependence in 
quantities on the dimensionless ratio m 2 / '-F 2 lat , as shown 
in Fig. [33J Similar analyses have been done before, for 
example, in Refs. [13,|9(J. (Note that f^.iat is replaced by 
/ir,iat(w/)/-\/2 at each m/ to be consistent with the chi- 
ral extrapolation formulation used in Refs. [87l. \8§L l89l|.) 
Since this is a dimensionless quantity, there is no system- 
atic error coming from the lattice scale determination; 
therefore, it is a good way to find out whether we have 
control over lattice artifacts. 

For the momentum fraction, (x) u —d, the fit yields 
an extrapolated value of 0.260(30) with an acceptable 
X 2 /dof = 0.48. That this extrapolation is consistent with 
the earlier ones indicates no further uncontrollable sys- 
tematics other than the extrapolation itself. The value is 
about three standard deviations above experiment. Note 
the larger error on the extrapolated point is caused by 
the "re-scaling" of the extrapolation range. Similar con- 
clusions are drawn for the helicity fraction, (x) Au-Ad- 
the updated extrapolation brings the value to 0.224(41), 
within two standard deviations of experiment. However, 
the x 2 per degree of freedom for this fit is still large at 
5.17. 

The fit results from the non-linear chiral extrapola- 
tions (see Appendix l 2 cl for more details) are summarized 
in Fig. [33J The momentum fraction is 0.171(19) at the 
physical point and agrees with experiment within one 
standard deviation. The fit has an acceptable \ 2 P er de- 
gree of freedom of 0.85. The helicity fraction extrapolates 
to 0.163(28) and is consistent both with the extrapolation 
in terms of m w and the experimental number. However 
the fit quality is poor: % 2 /dof = 4.57. These values are 
consistent with the earlier fits. 

The trial fits and extrapolations are summarized in 
Table [Xj In the cases of both the quark momentum and 
helicity fraction, the extrapolation with the chiral form 
seems to give us consistent values with the experiments. 
However, these extrapolations require a 50% drop from 
the magnitude of our measured points and do not al- 
ways give the smallest x 2 /dof. We simply must measure 
at lower pion masses before such extrapolations can be 
deemed trustworthy. Also note that the extrapolations 
in terms of m 2 or the dimensionless quantity m 2 /F% lat 
are consistent with each other. 

We emphasize that the data in general is not consistent 
with a simple linear form, but that it does trend toward 
experiment as the quark mass is reduced. 



D. Transversity 

Another interesting quantity regardin g th e spin struc- 
ture of the nucleon is its transversity [9 ll |92| . (l)s q . Here 
again we have calculated only the isovector quantity. The 
raw data arc shown in Fig. [34] together with the fit range 
we use; Table IXj summarizes the bare values. Again, at 
the mf = 0.02 point, we use two different source-sink sep- 
arations to control possible excited state contributions; 
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FIG. 33: Renormalized quark momentum fraction, (x) u -d, 
and helicity fraction, (iJab-aj, in terms of ml/ F% lat . The 
blue band indicates the linear extrapolation and the brown 
band is the phenomenological fit as described in Refs. [87], |88| . 
[8!|. (Graphics conventions as in Fig. [28]) 



in this case the two choices give results that are com- 
fortably consistent. Therefore, we use the value obtained 
from i scp = 10 for all further analysis and discussion. 

Using these fit values and the renormalization at the 
chiral limit listed in Table Hill we obtain the quark-mass 
dependence as shown in Fig. 1351 Similar to the axial 
charge radius, we observe a significant deviation at the 
lightest quark mass. We extrapolate linearly and obtain 
{VjSu-Sd = 0.93(6) + 0.37(17)m 2 with a x 2 per degree 
of freedom of 7.3. This translates into a value at the 
physical point of 0.93(6). It will be interesting to com- 
pare this result with experiment when the latter becomes 
available. 



E. Twist-3 moment 

We also calculated the twist-3 first moment of the po- 
larized structure function, d\. The good chiral symmetry 
of DWF prevents lower-dimensional operators from spoil- 
ing the calculation, unlike in calculations with fermion 
discretizations that violate chiral symmetry, such as Wil- 
son fermions. The calculated bare values are summa- 
rized in Fig. [35] an d Table Kill The quark-mass depen- 
dence of the bare values is plotted in Fig. [371 (Note 
that we do not intend to compare these calculations with 
experiment so do not renormalize them.) A linear fit 
yields -0.0059(16) + 0.126(4)m 2 and an extrapolated 
value di° — —0.0035(13) at the physical pion mass, 



which is about three standard deviations away from zero, 
but small compared to the values at non-zero mass. This 
suggests the Wandzura- Wilczek relation [6l| , which need 
not hold in a confining theory, nevertheless holds in QCD. 



IV. CONCLUSION AND OUTLOOK 

In this work, we have presented numerical lattice QCD 
calculations using two degenerate flavors of dynamical 
(DWF) quarks. We calculated isovector vector, axial, 
tensor and pseudoscalar form factors and some low-order 
moments of the unpolarized and polarized structure func- 
tions. 

We found that the ratio of the axial charge to vec- 
tor charge, gA/gvi has a significant drop at the light- 
est pion mass point (about 500 MeV) which may be a 
sign of finite- volume effect; the heavier pion mass points 
(around 600 and 700 MeV) are in rough agreement with 
experiment and show small dependence on the quark 
mass. A linear extrapolation in pion mass-squared gives 
0.89(6), which deviates from experiment by more than 
five standard deviations; leaving out the lightest point 
gives 1.23(12). If the downward trend at the lightest 
point, away from experiment, is indeed a finite- volume 
effect, then even m^L k 5 is not large enough for nu- 
cleon calculations with dynamical fermions, in contrast 
to quenched calculations [J] 2 . 



TABLE X: Summary of the extrapolations of the renormal- 
ized first moment of the unpolarized distribution. The chiral 
perturbation theory formulation can be found in Ref. [87t HH, 
|80| or appendix I 2 cl The fit parameters are in the order of 
fit slope, intersection for the linear fits and C, e for the chiral 
formulation. 





{x) u -d 


X 2 /d.o.f. 


fit parameters 


Linear vs. ml 


0.282(19) 


0.28 


{0.09(5), 0.280(19)} 


Linear vs. ml/F% lat 


0.260(30) 


0.48 


{0.0024(12), 0.26(3)} 


ChPT vs. ml 


0.147(9) 


2.19 


{0.109(8), 2.64(11)} 


ChPT vs. rr£/F% M 


0.171(19) 


0.85 


{0.139(17), 3.74(22)} 


experiment 


0.154(3) 








(a;) Au-Ad 


X 2 /d.o.f. 


fit parameters 


Linear vs. ml 


0.286(25) 


6.78 


{0.21(6), 0.281(27)} 


Linear vs. ml/F% M 


0.224(41) 


5.17 


{0.0061(17), 0.22(4)} 


ChPT vs. ml 


0.170(14) 


3.57 


{0.134(12),0.25(12)} 


ChPT vs. ml/Fl^ t 


0.163(28) 


4.57 


{ 0.135(27), 3.31(20)} 


experiment 


0.196(4) 







2 This effect may in fact be illusory due to the presence of quenched 
chiral logs as recently argued in |93H 
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TABLE XI: Bare transversity. 



m f 


tscp 


(1-}su-5d 


0.02 


12 


1.29(11) 


0.02 


12 


1.25(4) 


0.03 


10 


1.43(3) 


0.04 


10 


1.41(4) 



TABLE XII: Bare values of the twist-3 moment di of the 
polarized structure function. 



mf 


^sep 


ju—d 




4 


0.02 


12 


0.0240(25) 0.0185(19) 


-0.0055(13) 


0.02 


10 


0.0246(8) 


0.0189(7) 


-0.0056(4) 


0.03 


10 


0.0411(8) 


0.0325(6) 


-0.0085(3) 


0.04 


10 


0.0544(7) 


0.0428(6) 


-0.0117(3) 



We studied the momentum-transfer dependence of the 
Dirac, Pauli, axial and induced pesudoscalar form fac- 
tors. Using a conventional dipole extrapolation, we found 
masses My, and Ma at each pion mass points that 
are more than 50% higher than the experimental ones, 
while the mean-square charge radii, (r v ), (rf,), (r A ) 
are far below experiment. Similar effects have been ob- 
served in the past dynamical calculations Q • Most phe- 
nomenological studies expect to see a dramatic increase 
at lower pion mass; although this trend is not observed in 
our current study, it may be resolved in the future with 
lighter pion calculations. 

The magnetic moments of nucleon were calculated in 
two different ways: first, directly from extrapolation of 
the Pauli form factor and second, by looking at the 
ratio of the isovector electric and magnetic moments 
G v M (0)/G v E (0) in the forward limit; the latter expres- 
sion has milder momentum dependence. Using the ra- 
tio method, we see decreasing m 2 dependence, and with 
a simple linear extrapolation we find fi p — [i n — 3.4(7), 
roughly consistent with the experimental value. 

In studying the axial current form factors, we find that 
the renormalized axial form factor at finite momenta ap- 
pears to be in broad agreement with experimental val- 
ues from neutrino scattering. The Goldberger-Treiman 
relation appears to hold for low momentum transfer, 
q 2 < 0.5 GeV 2 . Assuming the relation, we obtain an es- 
timate for the pion-nucleon coupling of g^NN — 15.5(1.4) 
and g P = ^-G p (0.88to?,) = 7.68 ± 1.03, which is con- 
sistent with experiment. 

For the structure functions, the renormalized values 
of (x) u ^d and (x}Au-Ad show a trend toward their re- 
spective experimental values at the lightest quark mass 
point in our study. The quark momentum fraction, 
(^)«-d(2 GeV), is 0.282(19) by simple linear chiral ex- 



1.7 
1.5 
1.3 
1.1 

0.9 

1.7 
1.5 
1.3 
1.1 
0.9 



- ' 1 ' 








■ 


■ 






■ 


* 1 


1 1 






ni| 


=0.02, t 


scp — 


12 




I : 










» — c 








< 

> 




( 


) ( 


5— « 


>— « 




>— e 


> — C 








. . , . 






, 




, 








. , . j 



1.7 
1.5 
1.3 
1.1 
0.9 



6 

t-t, 



10 



12 



1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 

I m f =0.02, t sep =103E 

5 inxijiiifiis 1 




m f =0.04,t scp =10 



6 

t-t s , 



10 



12 



FIG. 34: The bare values of the transversity; symbols as in 
Fig. [27] 



trapolation, which overshoots the experimental value by 
more than three standard deviations. The quark helicity 
fraction, (x)Au-Ad(2 GeV), is 0.286(25) by simple linear 
extrapolation, which overshoots the experimental value 
by almost four standard deviations. The former does 
not deviate significantly from a linear ansatz when the 
source-sink separation is set to 10 time slices. When in- 
creased to 12, the lightest point lies below the straight 
line fit, suggesting possible excited state contamination 
for i scp = 10, and perhaps the desired physical quark 
mass dependence for £ sop = 12. However, the statistics 
for the larger separation are relatively poor, and defini- 
tive conclusions can not be drawn. A similar pattern 
holds for the helicity fraction, except that at t scp = 10 
a non-linear quark mass dependence is already observed. 
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FIG. 35: Nonperturbatively renormalized transversity with 
linear extrapolation to the physical pion mass. (Graphic con- 
ventions as in Fig. 1281 ) 



Non-linear fits motivated by chiral perturbation theory 
were also used to fit the data. Again, these indicated 
favorable trends in some cases, but definitive results re- 
quire simulations with lighter quark masses and multiple 
volumes, work that is now well underway. 

We also calculated transversity and twist-3 matrix 
elements. The transversity renormalized at 2 GeV, 
(l)j u -ta(2 GeV), is 0.93(6) by simple linear extrapola- 
tion in m^. However, there is noticeable non-linearity 
in the data. A linear extrapolation to the physical 
point yields di e = —0.0035(13) for the twist-3 operator, 
which is about three standard deviations away from zero, 
indicating only small breaking of the Wandzura-Wilczek 
relation 6 if . 



Unfortunately, even with our dynamical calculations, 
we cannot resolve the long-standing differences between 
the predictions of QCD and experimental measurements, 
in such quantities as the charge radii, and quark mo- 
mentum and helicity fractions. If both QCD and the 
relevant experiments are correct, we expect to see a dra- 
matic shift in these values as the constraints of high 
quark mass, small volume and coarse lattice spacing are 
lifted. The limitations of the current gauge ensembles 
do not allow such a study at this stage. However, we 
are currently studying these quantities with a new se- 
ries of ensembles produced by the RBC and UKQCD 
Collaborations [H, [86[ . And another ensemble is being 
generated by these groups and LHPC. They are respec- 
tively at a -1 = 1.7 and 2.1 GeV lattice cutoff with up 
to (2.75 fm) 3 volumes; they feature physical dynamical 
strange quark mass and two degenerate up and down 
quarks as light as 1/7 the strange mass. This will give 
us better control of the systematic errors in a calcula- 
tion with much lighter pion masses in the valence and 
the sea sectors; whether we will start to see the curva- 
ture suggested by the chiral perturbation theory or other 
phenomenological models will be very interesting to dis- 
cover. 
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FIG. 36: The bare twist-3 moment di; symbols as in Fig. 1271 
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APPENDICES 

1. Nonperturbative renormalization 

In this appendix, we describe the formulation we used 
in the nonperturbative renormalization in more detail. In 
subsection III CI we obtain renormalization constants in 
RI/MOM scheme. In order to compare our results with 
the experimental ones, we need to further to convert the 
renormalization constant to MS scheme. In the case of 
the first moments of the momentum fraction and helic- 
ity operators, one finds the conversion equations in the 
continuum via 168] 



'MS 



Z RI 



16tt 2 



G n + (1 - OSn-1 



4 2\ \E M P^(P) 



n+1 



with 



G n — 



n(n + 1) 
2 



(-3 - S„ 



25'„+i) 



(41) 



5,) 



1 

--4^-(2S r S H )-l, (42) 

j=2 J 
n 1 

(44) 



J 



3 = 1 

jU2,...,/l„ 

where £ = in Landau gauge. c, 

l 

3 



IMfJ,2...fl n 



■fay 



are 6^46^4 — | 53fc=i ^^kS v k for momentum fraction and 



<5/i3<5i/4 + 6^46^3 for helicity operators. In the tensor- 
current case, we use the matching factor calculated from 
Ref. M: 



z m 



= 1 



\4ttJ 81 



(-4866 + 1656C(3) + 259n/). 



(45) 



We use the three-loop solution for the running of 



Cts 
47T 



1 



/9i In In (/1 2 /A 



QCD J 



/3 ln(/i 2 /A 2 QCD ) /3 3 ln 2 (>/A 



2 

QCD 



/3 5 ln 3 (m 2 /A 2 



QCD 



ft 2 In In (mVA^cd) 



{/3 2 ln 2 ln ( M 2 /A 2 QCD ) 



(46) 



where a = |- and Cf = 
The various & are 



00 
347V 2 



2N C 
llNr - 



A QCD = 300 MeV[7C 
2n f 



2857 
~54 
1415 , 

11 TV 2 - 1 



4/V 2 



79 



-n f 



N, 

205 
~36~ 



(#c ~ l)n/ 



—N c ni in MS scheme, (47) 



-r Jtx -r j 

18 iV c 7 54 

where 7V C is the number of colors and n/ is the number of 
quark flavors. In the quenched approximation, rif = 0; 
in this work, 71/ = 2. 

The running of the Z-factors to two loops [§5| is 

An 



7i - Mo) , (48) 



where 7^ is defined by the anomalous dimensions 7$ di- 
vided by 2(3q. The anomalous dimensions for twist- 
two operators of interest were calculated by Floratos 
et al. 691: 



7o 



7i 



2C f 



n(n + 



(-8.30C 2 , + 27.85C F C A - 9.48C F 7», (49) 



for twist-n operator. (Note that where Tp(n,f) 



2N, 



- and Ca = N c .) The anomalous dimensions 



C F = 

for the operator, such as O34 , can be found in Ref. [96 



70 = 2C F 
Cf, 



7i 



-£-(257C A - 171Cir - 522». 



(50) 
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2. Chiral extrapolation models 

In this appendix, we collect some relevant formulae 
used in this paper from various models proposed to con- 
duct chiral extrapolations for baryonic observables. 



a. Finite-volume correction on qa 

The small-scale expansion scheme [72[ tries to use ex- 
plicit degrees of freedom from the pion, nucleon and 
A(1232) resonance, expanding in terms of A , the mass 
splitting between the N and A in the chiral limit. This 
splitting is treated as a small parameter, of the order 



0(e). They define g A and F° as the value of the axial 
charge gA and the pion decay constant in the chiral 
limit, ca the iVA and gi the A A axial coupling constants 
respectively. We simply abbreviate A , g A and as A, 
gA and F® in the following formula. 

At the leading order 0(e 3 ), the finite- volume correction 
to the nucleon axial charge, gA, was proposed in Ref. (7lj | 
within the SSE scheme: 



9aF 3 H 



9aFi 
2c 2 A g A F 4 



2ci 
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9 A - ^.91 



(51) 



where 



Fi(m,L) 

F 2 (m,A,L) 
F 3 (m,L) 

F 4 (m,A,L) 



E 

n=£0 



K (mL\n\) - 
Ki(mL\n\ 



Ki(mL\n\) 
mL\n\ 



min 



E 

3 ^ Xi(m£|n|) 
2 V mL\n\ : 

n^O 1 1 



-K {mL\n\) 



A 



dfi 



2l3K (/3L\n\) + (A 2 - m^ln^^nl) 
y/(3 2 + A 2 - m 2 



E 

n#0 L 



Ki{mL\n\) ire~ mL ^ A 2 - 



mL\n\ 



2AL|n| 



>A 



d/3 



13 K (J3L\n\) 
v//? 2 + A 2 - m 2 



(52) 



with modified Bessel function of the second kind K a {z). 
(Note that the m in the above equations are short 
for m^.) This formula can be reduced to that of LO 
HBChPT with c A = 0. 

b. The mean-squared Dirac radius in HBChPT 

Heavy baryon chiral perturbation theory 
(HBChPT) 97] involves only the pion and nucleon, 
expanding in terms of the momentum p. At leading 
order (LO), once we use the physical values of gA, F„ 
and {r v ) instead of their values in the chiral limit, the 
HBChPT formula provides a closed form for the m 2 
dependence of the Dirac mean-squared radius without 



unknown parameters [zD, [9?| : 

1 -|- §g A ( \ 

(r v )(m« M ) = (r v ) exp - J^r^^ ) (53) 

where the scale /i should be set by the physical value of 
the pion mass as fi = 0.139 GeV 3 . 



c. Chen, et al. 

Chen et al. [13, [H, 89] calculate a chiral perturbation 
expression in the continuum for the moments of the un- 
polarized and helicity distributions: 
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3 A missing factor of 2 in the prefactor of In rra^ term, which can be quark mass region, 

found in Ref. 1981 . causes an accidental agreement between this 
leading one-loop form and lattice data obtained in the heavy 
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(54) 
(55) 



We use /x = 0.139 MeV and F v , e ^ = 92.4 MeV in our 
extrapolation. 

There is a variation of this formulation replacing the 
^7r,exp with F^iat and changing the scale to i at . This 



provides a dimensionless quantity F £' lat in the extrapo- 

lation, which is independent of the scale setting in the 
lattice calculation. 
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